Quantum Hall plateau transitions in disordered superconductors 
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We study a delocalization transition for non-interacting quasiparticles moving in two dimensions, 
which belongs to a new symmetry class. This symmetry class can be realised in a dirty, gapless 
superconductor in which time reversal symmetry for orbital motion is broken, but spin rotation 
symmetry is intact. We find a direct transition between two insulating phases with quantized Hall 
conductances of zero and two for the conserved quasiparticles. The energy of quasiparticles acts as 
a relevant, symmetry-breaking field at the critical point, which splits the direct transition into two 
conventional plateau transitions. 



The variety of universality classes possible in single- 
particle models of disordered conductors is now appreci- 
ated to be quite rich. Three of these classes were identi- 
fied early in the development of weak localization theory 
PJ : they are distinguished by the behavior of the system 
under time-reversal, and by its spin properties, and are 
termed orthogonal, unitary and symplectic, in analogy 
with Dyson's classification of random matrice ensembles. 
Further alternatives can arise by two different mecha- 
nisms. First, in certain contexts, most notably the inte- 
ger quantum Hall effect, the non- linear a model describ- 
ing a two-dimensional system may admit a topological 
term Q , which results in the existence of extended states 
at isolated energies in an otherwise localized spectrum. 
Physically, such systems have more than one distinct in- 
sulating phase, each characterised by its number of edge 
states, and separated from other phases by delocalization 
transitions. Second, it may happen that the Hamilto- 
nian has an additional, discrete symmetry, absent from 
Dyson's classification. This is the case in two-sublattice 
models for localization, if the Hamiltonian has no matrix 
elements connecting states that belong to the same sub- 
lattice HQ . It is also true of the Bogoliubov-de Gennes 
formalism for quasiparticles in a superconductor with dis- 
order [^|^] . One consequence of this extra symmetry is 
that, at a delocalization transition, critical behavior can 
appear not only in two-particle properties such as the 
conductivity, but also in single-particle quantities, such 
as the density of states. 

Universality classes in systems with extra discrete sym- 
metries of this kind have attracted considerable atten- 
tion from various directions. A general classification, 
systematizing earlier discussions JIJ5), has been set out 
by Altland and Zirnbauer ||, who examined mesoscopic 
normal-superconducting systems as zero-dimensional re- 
alisations of some examples. Very recently, quasiparticle 
transport and weak localization has been studied in dis- 



ordered, gapless superconductors in higher dimensions, 
with applications to normal- metal/superconductor junc- 
tions and to thermal and spin conductivity in high 
temperature superconductors p[-^0|. Separately, the be- 
haviour of massless Dirac fermions in two space dimen- 
sions, scattered by particle-hole symmetric disorder in 
the form of a random vector potential, has been inves- 
tigated intensively as a tractable example of a dis- 
ordered critical point. And much before this, the one- 
dimensional tight binding model with random nearest- 
neighbour hopping was shown jl^] to have a delocaliza- 
tion transition and divergent density of states at the band 
center, the energy invariant under the sublattice symme- 
try. 

In this paper, we study a new delocalization transi- 
tion in two dimensions that combines both of the above 
features: the transition separates phases with different 
quantized Hall conductances for the quasiparticles, and it 
occurs in a system which has a discrete microscopic sym- 
metry. This transition can take place in a gapless super- 
conductor under appropriate conditions: time-reversal 
invariance for orbital motion must be broken by an ap- 
plied magnetic field, but the Zeeman coupling should be 
negligible, so that the full spin-rotation invariance re- 
mains intact. Formally, we suggest that the model whose 
behavior we examine numerically is a representative of 
the symmetry class labeled C by Altland and Zirnbauer 
, and that the delocalization transition is associated 
with a topological term allowed in the field theory, as 
noted by Senthil et al The possibility of quantum 
Hall states in superconductors with broken parity and 
time-reversal symmetry has been emphasised by Laugh- 
lin [ ft3| . A direct transition into such a phase, in the 
presence of disorder, is of particular interest in connec- 
tion with the theory of the quantum Hall effect, since 
it is between phases with Hall conductance differing by 
two units. Changes of Hall conductance by more than 
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one unit at a derealization transition are precluded in 
generic systems, by the standard scaling flow diagram 
fuj for the integer quantum Hall effect, and are possible 
in the model we study only because of additional sym- 
metry. In the presence of a symmetry-breaking coupling 
in the model, of strength A, the transition occurs in two 
stages with a separation that varies as A* for small A, 
where (f> ~ 1.3. Such a coupling is introduced if one treats 
quasiparticle motion at finite energy. 

The system we consider is formulated as a generali- 
sation of the network model jl5| for the quantum Hall 
plateau transition. The original version of this model de- 
scribes guiding center motion of spin-polarised electrons 
within one Landau level of a disordered, two-dimensional 
system in a magnetic field. It therefore has broken time 
reversal symmetry for orbital motion, and contains no 
spin degree of freedom. It is specified in terms of scat- 
tering or transfer matrices, defined on links and nodes 
of a lattice. The model can be generalised in various 
ways. Spin can be incorporated by allowing two ampli- 
tudes to propagate on each link. This has been done pre- 
viously |l6|-|l8|], with the intention of describing a spin- 
degenerate Landau level in which the two spin states 
are coupled by spin-orbit scattering. In that case, the 
random U(l) phases which characterise propagation on 
the links of the original model are replaced with ran- 
dom U (2) matrices, mixing the two spin states without 
any rotational symmetry. In the work presented here, we 
choose instead random SU(2) matrices, preserving spin 
rotational symmetry. 

In detail, the transfer matrix associated with each link 
of the model is an SU(2) spin rotation matrix of the form 
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where 5i,S2,x are random. The transfer matrix at the 
nodes is parameterized by e ± ^A so that the transmis- 
sion probability for the two spin states is [l + exp(— 7r(e± 
^A)] -1 respectively. The value of e determines the Hall 
conductance of the system, as measured at short dis- 
tances: varying e drives the model through the dereal- 
ization transition. A non-zero value for A breaks spin- 
rotation invariance, and will in fact change the universal- 
ity class for the transition. Collecting factors, the trans- 
fer matrix across one node and the links connected to it 
is a 4 x 4 matrix of the form fllCI 
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and Ui are as given in Eq (Q) . From these 4x4 transfer 
matrices, T, one can build up a larger transfer matrix, 
of size 2Mi x 2M;, with Mi even, to describe scattering 
in one slice of a system of width Mi links, (which has 
M = 2Mi scattering channels) by using independent re- 
alisations of T as diagonal blocks of the larger matrix. 

Both the 4x4 transfer matrix, T, and the larger ones 
derived from it, which we denote here also by T, are 
invariant under an antiunitary symmetry operation rep- 
resenting spin reversal. The corresponding operator is 
Q = TL<E>iT y ■K 1 where the Pauli matrix r y acts on the two 
spin states propagating along each link, 1 is the Mi x Mi 
unit matrix, and K is complex conjugation. The symme- 
try, QTQ 1 = T, holds for A = only; it implies that 
the Lyapunov exponents of the transfer matrix product 
have a two-fold degeneracy at A = 0, which we exploit 
in the analysis of our simulations, as described below. 

It is possible to relate this network model to a Hamil- 
tonian, H, following Ref [lSf! , by constructing a unitary 
matrix which can be interpreted as the evolution oper- 
ator for the system, for a unit time-step. Taking the 
continuum limit, one obtains, in the case of the original 
network model for a spin-polarised Landau level, a two- 
component Dirac Hamiltonian with random mass, scalar 
potential and vector potential. For the network model of 
current interest, we get instead a four-component Dirac 
Hamiltonian of the form 



H = (<r x p x + (J z p z + may) (8 1 + 1 ® ol 



(4) 



where <Ji and Ti for i = x,y, z are two copies of the Pauli 
matrices, 1 is the 2x2 unit matrix, p x and p z are the two 
components of the momentum operator in the plane of 
the system, the mass m is proportional to e, the distance 
from the critical point, and the real, three-component 
vector, a, is a random function of position. This Hamil- 
tonian has the symmetry QHQ 1 = —H p0| , which is 
the defining feature of the class labelled C by Altland 
and Zirnbauer || . A non-zero value for A in the network 
model introduces an additional term, H ' = A a y ®T z into 
the Dirac Hamiltonian, breaking the symmetry. Equally, 
since the symmetry relates eigenstates with energies ±E, 
and leaves invariant only those at energy E = 0, non-zero 
E, like A, acts as a symmetry-breaking perturbation. 

The models represented by Eqs @ and (^) describe 
propagation of quasiparticles which are conserved, and 
which are obtained within the Bogoluibov-de Genncs 
formalism by making a particle-hole transformation on 
states with one spin orientation (see, for example, Ref 
§ ) . Specifically, starting from the Bogoluibov-de Gennes 
Hamiltonian for a singlet superconductor 
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and introducing transformed operators, 7^ = Cj-f and 
7], = , one has 




This Hamiltonian, like the SU(2) network model that we 
simulate, has the symmetry QHsQ 1 = —H,s [0], where 
here Q = ir v K and r y acts on the particle-hole spinor of 
Eq. (5). For singlet pairing (maintaining spin-rotation 
invariance) Ay is symmetric while hij is hermitian and 
the symmetry under Q is obvious; time-reversal symme- 
try is broken if h 7^ h* or A ^ A*. Since quasiparticle 
charge in a superconductor is not conserved, and charge 
response is anyway controlled by the condensate, the lo- 
calization problem of interest in a disordered supercon- 
ductor, as emphasised in Ref ||, involves spin and energy 
transport, rather than charge transport. We stress that 
the Hall conductance examined below is a property of 
quasiparticles described by a Hamiltonian such as Eq (||). 

We study the model defined from Eq (§) at a range of 
values for e and A. Preliminary calculations, reported 
earlier pH , were limited to A = 0. We compute the nor- 
malized localization length, £,m/Mi, for strips of width 
AI = 16, 32, 64, with periodic boundary conditions. For 
small A we find it necessary also to use M = 128, 256, 
in order to identify more clearly the critical properties. 
The parameters Si/2tt, ^2/277, x are each chosen to be uni- 
formly distributed in the interval [0, 1], so that each ma- 
trix U associated with a link is independently distributed 
with the Haar measure on SU(2). Runs were carried out 
for strips of length 60,000 (for M = 16, 32), 240,000 (for 
M = 64), and 480,000 (M = 128,256). The errors are 
typically 0.5%, except for M — 256. 

The behavior at non-zero A is shown in Fig. 1, region I. 
For A = 2.0, extended states (£,m/Mi independent of M) 
appear clearly at two energies, ±e c , with e c (A = 2) w 0.6. 
A one-parameter scaling fit for £,m/M = /[(e— e c )M 1 / I/0 ] 
yields isq w 2.5, the conventional quantum Hall exponent 

0,011. For A = °- 2 ' £ c is t0 ° sma11 t0 be resolved 
by this method (Fig. 1, region II). Proceeding in this way 
at a range of values for A, we construct the phase dia- 
gram for the model, shown in Fig. 2. With A / 0, the 
two phases at e < —1 and e ^> +1 are separated by 
an intermediate, small e, phase. Counting edge states 
in each phase in the strong-coupling limit (e — > ±00 at 
fixed A, and A — > 00 at e = 0, respectively), we find 
that the quasiparticle Hall conductance takes the values 
0, 1, and 2 in successive phases with increasing e. If A is 
made smaller, the boundaries of the intermediate phase 
approach each other, and at A = there appears to be a 
direct transition between phases with Hall conductance 
differing by two units. 

To study this direct transition, we examine behavior 
at small A in more detail. On the line A = 0, the local- 
ization length diverges at a single critical point, e c = 0, 



with an exponent, v = 1.12, which is different from that 
at the conventional plateau transition (Fig. 3, curve I). 
Close to this line, scaling with system size is quite com- 
plex and, in particular, the variation of £m/M with M 
is not monotonic. In order to extract scaling proper- 
ties at small, non-zero A, we monitor the deviation from 
Kramer's degeneracy of the smallest two positive Lya- 
punov exponents, Ai and A2, of the transfer matrix prod- 
uct that represents the sample, defining £ = M(X 2 — Ai). 
Finite size scaling of both £ and £m is shown in Fig. 3 
(curves II and III), as a function of A along the symme- 
try line e = 0. We find for deviations from Kramers de- 
generacy £ = /i(AM, ) and for the localization length 
£m/M = / 2 (AM, 1/Al ), with fj, w 1.45. The maximum in 
/2 indicates that states at e = 0, A 7^ appear delocal- 
ized at small M while at larger M the true, large-distance 
localized nature becomes apparent. Thus to see localiza- 
tion one needs M > (9/A) 145 : for example, for A = 0.2 
one needs M > 200. 

We propose, then, that e = A = is a critical point 
at which e parameterizes the symmetry-preserving rele- 
vant direction, and A is a symmetry-breaking field, so 
that !;m/Mi is described near the fixed point by a two 
parameter scaling function 

^ M /M l = f(eMl /u ,AM^). (6) 

with v = 1.12 and /1 as 1.45. In the presence of a 
symmetry-breaking field, A 7^ 0, scaling flow is away 
from the new fixed point, giving quantum Hall plateau 
phases except on trajectories which connect this unstable 
fixed point to fixed points at finite A, representing the 
conventional universality class for plateau transitions. At 
these, a finite critical e c (A) is expected with an exponent 
v , so that £ ~ [e - e c (A)]- u ° . Since the values of eM 1 /" 
and AAf 1 / M on a critical trajectory serve to define a one- 
parameter curve, we expect e c (A) = zLcA^^ . Thus, as 
A approaches zero, extended states coalesce, having a 
separation, 2e c cx A 13 (the line in Fig. 2), which is much 
smaller than A, their separation in the absence of cou- 
pling between the two spin orientations. 

A further aspect of the critical point which is of inter- 
est, but not accessible within our numerical approach, is 
the behavior of single-particle quantities such as the den- 
sity of states, discussed recently in Ref [10]. We expect 
for the Hamiltonian of Eq ([|) a finite density of states at 
all energies provided A/0, and singularities in the den- 
sity of states at zero energy when A = 0, with a different 
nature according to whether e = or e 7^ 0. 

In conclusion, we have shown that quantum Hall 
plateau transitions belonging to a new universality class 
occur in a model for a gapless superconductor which is in- 
variant under spin rotations, but which has time-reversal 
symmetry broken for orbital motion. In contrast to the 
conventional plateau transition, the Hall conductance for 
conserved quasiparticles changes at this transition by two 
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units. We have examined critical behavior, and shown 
that there is a symmetry-breaking perturbation which is 
relevant at the critical point, splitting the transition into 
two, with extended states that coalesce as the symmetry- 
breaking field is removed. 
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FIG. 1. Normalized localization length £m/Mi. The symbols O, +, □, x, and A correspond to system widths M of 16, 32, 
64, 128, and 256 respectively. The symmetry breaking parameter is A = 2 (region I) or A = 0.2 (region II). 
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FIG. 2. Phase diagram. The + symbols indicate fitted positions of extended states, and the lines are e = ±0.25A 1,3 . 
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FIG. 3. Scaling functions: (I) Normalized localization length £m/Mi as function of eM 1 '" with v = 1.12 for A = 0. (II) 
Deviation from Kramer's degeneracy £ as function of AM 1 '" with // = 1.45 for e = 0. (Ill) £ M /M; as function of AM 1 '" with 
/i = 1.45 for e = 0. Symbols denote system widths as in Fig. 1. 
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